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We investigate the interaction of excitons in a two dimensional lattice and photons in a planar 
cavity in the presence of disorder. The strong exciton-photon coupling is described in terms of 
polariton quasi-particles, which are scattered by a disorder potential. We consider three kinds of 
disorder: (a) inhomogeneous exciton energy, (b) inhomogeneous exciton-photon coupling and (c) 
deviations from an ideal lattice. These three types of disorder are characteristic of different physical 
systems, and their separate analysis gives insight on the competition between randomness and light- 
matter coupling. We consider conventional planar polariton structures (with excitons resonant to 
photon modes emitting normal to the cavity) and Bragg polariton structures, in which excitons in a 
lattice are resonant with photon modes at a finite angle satisfying the Bragg condition. We calculate 
the absorption spectra in the normal direction and at the Bragg angle by a direct diagonalization 
of the exciton-photon Hamiltonian. We found that in some cases weak disorder increases the light- 
matter coupling and leads to a larger polariton splitting. Moreover, we found that the coupling of 
excitons and photons is less sensitive to disorder of type (b) and (c). This suggests that polaritonic 
structures realized with impurities in a semiconductor or with atoms in optical lattices are good 
candidate for the observation of some of the Bragg polariton features. 

PACS numbers: 71.36.+C, 78.20.Bh, 78.67.Hc, 71.55.-i, 78.40.Pg 



I. INTRODUCTION 

There is currently a considerable interest in the physics 
of strongly coupled light-matter systems. Examples in- 
clude experiments on the Bosc-Einstcin condensation 
(BEC) of microcavity polaritons ^L^A which are mixed 
states of excitons and photons (see Refs. [H and @) 
and on the strong coupling of cavity photons and quan- 
tum dot excitons In order to combine the prop- 
erties of matter states in a lattice and planar photon 
modes, we have recently investigated a structure con- 
sisting of an array of quantum dots in a planar cav- 
ity^ as schematically shown in Fig. [TJ This particu- 
lar geometry for light-matter coupling can be realized 
in many different ways using either semiconductor-based 
systems (quantum dots, impurities, metallic gates) or 
atomic systems, such as atoms in optical lattices. Exper- 
imental and theoretical investigations show that impu- 
rity bound excitons can have very small inhomogeneous 
broadening and strong confinement energy (e.g. of the 
orders of 50 me V 11 ' 12 ). Moreover, ordered arrays of sin- 
gle dopants,— as well as the control of a single impurity 
using scanning tunneling microscopy have been experi- 
mentally demonstrated. 14 i 15 

Polaritons in structures similar to the one in Fig. [TJ but 
involving photonic crystals^ or optical lattice a 17 ' 18 have 
been recently investigated theoretically. Many of these 
structures have planar spatial periodicity, which leads to 
very interesting and novel properties. We have recently 
explored theoretically Bragg polariton modes at some 
special symmetry points of the Brillouin zone bound- 
aries^ These zone-edge Bragg polaritons can have ex- 
tremely small effective masses, typically three orders of 
magnitude smaller than conventional cavity polaritons, 
and behave effectively as Dirac quasiparticles, similar to 




FIG. 1: Scheme of the system studied. A lattice of quan- 
tum dots or impurities is embedded in a planar microcavity. 
The energy, position, and oscillator strength of the localized 
excitons fluctuates from site to site. 

light-mass (relativistic) electrons in graphene, which have 
been experimentally investigated recently 19 ' 20 Polaritons 
with Dirac dispersion could have interesting applications 
to the physics of polariton BEC and superfluidity, as well 
as for spin-coupling control^ 

In Fig. [TJ we show disorder effects due to fluctuations 
in the quantum dot size and position. These disorder 
effects are the main focus of this paper. The effect of dis- 
order on polaritons has been extensively studied in the 
quantum well-microcavity case (see Ref. |5j and references 
therein). It was found that if the potential fluctuations 
due to the disorder are comparable with the Rabi split- 
ting then the two polaritonic peaks disappear and one 
inhomogeneously broadened peak remains in the absorp- 
tion spectrum. Furthermore, a similar behavior has been 
theoretically predicted for an ensemble of two-level sys- 
tems coupled to a single photon mode^ 2 - The dynamics 
of the wave packet and its localization in disordered one- 
dimensional cavity has been been studied too.— In this 
paper, we investigate the effect on polariton modes of en- 
ergy, oscillator strength, and position fluctuations. The 
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numerical results are discussed and compared to a per- 
turbation theory approach. The paper is organized as 
follows: The theory is introduced in Sec. II. followed by 
results and discussion in Sec. III. Conclusions are drawn 
in Sec. IV. 



II. THEORY 

We are going to investigate a disordered planar lattice 
embedded in a planar microcavity structure. The role of 
the lattice is to localize excitons and the physical mecha- 
nism of localization will not be specified at this point. As 
discussed above, quantum dots, impurities, metal gates 
and optical lattices can be used to localized the exciton 
states. We further assume that only one excitonic level is 
present on each site, implying a strong localization. This 
also allows us to describe the exciton transfer process be- 
tween sites using a tight-binding approach. The case of 
quantum dots has been extensively discussed in Ref. [ToL 
We start with the Hamiltonian in second quantization 
(h = 1) 

h = -,(■:(■; • /x E c) +k c 3 

j j.keNB 
9 31 

where Cj is the exciton creation operator on the j th site 
at position Rj with energy ujj and exciton-photon cou- 
pling gj q , tx is the inter-site energy transfer between 
nearest neighbors sites, a q is the creation operator of the 
cavity photon mode with in-plane momentum q and en- 
ergy u) q . By writing the coupling constant as gj q e lqRj 
we can analyze separately the disorder effect induced by 
position fluctuations and oscillator strength fluctuations. 



A. Disorder properties 

We assume that the distance between localization sites 
is much larger than the effective exciton localization 
length on each site. This justifies why fluctuations are 
site-uncorrelated. For energy fluctuations, defined as 
Aujj = ujj — u>x, ^->A" being the average exciton energy, 
we have 



and 



(Aojj) = 



(Aw m Aw„) = a^Sr, 



(2) 



(3) 



where (.) indicates averaging over a Gaussian ensemble. 
Similarly, we have for the oscillator strength disorder 

(9jq) = (dog) = 90q~, {9kq9lq') ~ {9kq){9lq') = o\ q ^kl ,(4) 



where we assume go q — ge q 13 / 4 with the characteristic 
size of the site^S (3 and 



a 2 , - a 2 e- ( - q2+q ' 2)02 / i 
°qq ~ g c > 



(5) 



where a 2 = (g 2 ) — (g) 2 is the on-site oscillator strength 
variance at q = 0. Finally, we obtain for the positional 
disorder 



(Rj) = Rj, (R k Ri) - (R k )(Ri) = vlku 



/ e ig-Rj\ _ e iqRj & -q ct^/2 



(6) 



where the positions Rj identify the ideal two dimensional 
lattice. Furthermore, we assume for simplicity that the 
different kinds of disorder are uncorrelated. This allows 
us to treat them separately in the perturbation approach 
described below. 



B. Polariton scattering 

We can separate the fluctuation-independent part of 
the Hamiltonian by defining the Fourier transform of the 
exciton operators as 



Cj = -L y e-^c, 

V qGl.BZ 



1 



q> 



(7) 



where Rj corresponds to the j site of a two dimensional 
ideal lattice and N is the total number of sites. The 
sum over q states is restricted to the first Brillouin Zone 
(BZ) of the reciprocal lattice space. This fluctuation- 
independent Hamiltonian then takes the form 

H = y^(u>xgClC q + u q+ Qa\ +Q a q+ Q 
q Q 

+ Y,^q+Q a l+ Q Cq + h.c.)Y (8) 



where Q is a reciprocal lattice vector, g q = VNgo q is the 
renormalized coupling constant, and the exciton disper- 
sion of a square lattice reads 

ux q = u x - 2t x (cosfeo) + co8(q y a)) , (9) 

where a is the lattice constant. We can then define the 
disorder coupling constants 

Au>„ = —=y^e iqR 'Aco4, 

q \/N ^ 

j 

= -^E ei(W) ^-5o 9 ), 



E( 



MR j)e -iq'R S g^ (10) 
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where the index o (p) indicates oscillator strength (posi- 
tion) disorder. The statistical properties of these func- 
tions are determined by their average value as 

(Aw*) = 0, 

(»&) = o, 

lqq> 



K 



(11) 



and correlations according to 

(Vqq'Vkk') = ^2&q-k,q>-k>+KOqk, 



(12) 



K 



Kq'^kk') ~ (Vqq>)(Vk%>) =Yj S 4- k >4'- k '+ K ^ k > 

K 

where K is a reciprocal lattice vector and 
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1) 



(13) 



£q = 90 q (e 

C qk = 5o 9 55 fe ( e - k+fe|2 ^ /2 -e- 92 ^/ 2 e- fc2 ^/ 2 ). 

The disorder terms for the three different mechanisms 
can then be written in a compact form as 

V qk 



h 0{p) = £K?4^' + ^-). 

qq' 



(14) 



Ho\P n q) — A-nq\Pnq) 



We can find the eigenvalues and eigenvectors of the dis- 
order free Hamiltonian by solving 

(15) 

These are the disorder-free polariton states that can be 
written in the form 

\Pnq) = Plq\0) = (Un q Cl + Y,V nq+Q al +Q )\0),(l6) 

Q 

where n is a band index, v nq and u nq are Hopfield 
coefficients^ and |0) is the exciton-photon vacuum. We 
can use these states and obtain an effective disorder po- 
tential for polaritons as 

Vnn'qq' = {P<nq\H l(2)\Pn' q>) ■, 

where J € {e, o,p} labels energy, oscillator strength, and 
positional disorder, respectively. This gives for the en- 
ergy disorder 



V e , , 

v Tin' qq 



N 



(17) 



and for the oscillator strength (positional) disorder 

yo(j>) . V"* f„°(p) „ , ,„* i *°(P) ..* ,,* ^ 

v nriqq' " 2—1 \ 'l+Qg' n 1 U nq+Q ^ 'Iq+Qq' u nq u n' q' +Q J 



(18) 



C. Absorption spectrum 



The absorption spectrum at an angle determined by 
the inplane q of the cavity photon can be calculated using 



the full propagator of the disordered system, which can 
be written in terms of the disorder-free polariton states 
as 



nn'kk' 



(19) 



with 



G n n'kk'{u) = G nk (L))5 nn '5kk' 

+ G° nk (co)T nn , kk ,(u)G° n , k ,(u), (20) 



where 



1 



oj - A nk + ie 



(21) 



is the Green's function for disorder-free polaritons. T is 
the scattering T-matrix that can be expressed as 

T(w) = ^T«( W ), 
j 

T^LM = (Pnk\V(G°(u)V) j \Pn'k>). (22) 

The imaginary part of the polariton propagator projected 
on a photon mode at a given wavevector I7) = |fc 7 ) gives 
the absorption spectrum at the corresponding emission 
angle as 



A 1 (uj) = -ImG 7 (w) , 



(23) 



with 



G 7 M = (7|G(w)|7) = £^nfe T <'fe 7 G nn r k ^{w). 



We can explicitely average over the disorder configura- 
tions to obtain 



with 

(Gwk^M) = G° nk ((j)5 nn ,5 kk , (24) 
+ G°,>)(T„„^ fc >))G°, fc >). 

In the case of disorder due to energy inhomogeneity, the 
ensemble averaging leads to the conservation of the mo- 
mentum, {T nn > kk /{w)) = {T nn > kk (uj))5hk' , since the aver- 
aging selects only diagonal terms in Eq. (fT2")h We note 
that this is the same as in the case of one particle scat- 
tering on an impurity ensemble.— The diagonal elements 
of the T-matrix can be interpreted as an effective self- 
energy, T, nk (u>) = T nnkk (u)) and within the pole approxi- 
mation H n k(to = A„fc), their real part gives a shift of the 
polariton levels. The shift can also be calculated using 
a direct perturbation theory approach as shown in the 
following section. 



J 
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D. Perturbation theory 

Here we calculate the corrections to the energy up to the second order in the disordered potential, using the 
statistical properties of the disorder derived above. We start with the first order energy shift for all three kinds of 
disorder 



E^ = ^\u nk \ 2 , E^ = 2Re {e%?W<* + * \ , (25) 



'N 

where u n k and v* n k+K are defined in Eq. ([14")) . The average value of the energy shift 



= 0, (26) 
(E'J) = 0, (27) 

(28) 



( E nk) = 2Re |^^+ KW " fcU ™fc+^| ' 



and variances are 



mm = -fM-, (29) 

({El'D 2 ) = 4Rc | E ^+K. k+ K' \u nk \ 2 v; ik+K v* nk+K , I , (30) 
((E'J) 2 ) (E^) 2 = 4Rc | Ck + Kk + K'\u nk \ 2 v* nk+K v* nk+K , \ . (31) 



K,K' 

The second order energy shift in the non- degenerate case is given by 



n k 

E ( o(p) * *o(p) * * \|2 

K\nk+Kk ,U rik' v nk+K + Vk+Kk' U nk V n>k'+K)\ /„ Q \ 

^nk ~ A A " ^ 66 > 

n' k' 

In the case of energy disorder, average and variance can be calculated as 

- \Unk\ 2 \u n 'k'\ 2 



( E nk) = -ftfL A_,._ A_„, ' (34) 



N ^ Ank ~ An'k' 
n k 

//Ei2,e-,2\ / 7^2, ex2 a t I U nk | 4 \u n 'k' | 4 , Q _v 



The average value for the the oscillator strength and positional disorder are given by 

°(p) 

( E nk iP) ) = E A nk l n T ■ ( 36 ) 
, , , **-nk ^n' k' 
n' k' 

where 

KK' 

+ K-kKK' U nk u *n'k' v nk+K'V n 'k' + K + $1°^ K ,\u nk \ 2 V n , k , +K V n , k , +K ,^ , (37) 

,q 1o (p) . (°{V) *o(p) K 

u kk'KK< ~ \'lk+Kk''lk+K'k'h 

q2o(p) . ; o(p) o(p) , 

u kk'KK' — \'lk+Kk''lk'+K'k/- 

Expressions for the variance are not shown but their calculation is straigthforward. The degenerate case for the second 
order shift will be discussed in the next section. ^ 

III. RESULTS AND DISCUSSION the case of an excitation at the M point (finite excita- 



Here we present the absorption spectra calculated nu- 
merically for various strengths of disorder. We consider 
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FIG. 2: (Color online) Polariton absorption spectrum in log- 
arithmic scale for (i) energy disorder with variance from 
(black) to 5 meV (magenta) for polaritons at two high sym- 
metry points of the first Brillouin Zone: M (a) and V (b). A 
homogeneous Lorentzian broadening of 7=50^eV was consid- 
ered. Material parameters correspond to ZnSe/CdSe systems, 
and can be found in Ref. 03. 

tion angle corresponding to qm = ( 7T /o-i' K / a )) an d at the 
r point (normal excitation). The exciton energy is tuned 
into resonance with the corresponding cavity modes. We 
have used material parameters of CdSe/ZnSe systems as 
summarized in Table III of Ref. Oil The full Hamilto- 
nian corresponding to a lattice with 31 2 sites and photon 
modes in an energy window [0.9u>Xi l-la>jf] is diagonal- 
ized. In order to make the results independent of the 
system size we keep fixed the total disorder-free coupling 
constant g q and we use periodic boundary conditions. 

Let us start with the effect of the energy disorder in 
the absorption spectra as shown in Fig. O The disorder- 
free case (black line) shows lower (LP) and upper (UP) 
polariton peaks with half-exciton and half-photon char- 
acter. Additionally, in the M -point case (a), there is a 
central purely photonic peak. The polariton splitting, i.e. 
the energy difference between the lower and the upper po- 
lariton energy, is larger by ~ 0.5 meV at the T-point than 
at the M-point due to the stronger coupling in the nor- 
mal direction^ The disorder shifts the LP (UP) energy: 
According to the perturbation theory the first order shift 
of the LP (UP) is zero (as seen from Eq. (|26p ) and the 
second order shift is negative (positive), according to Eq. 
P4p . and leads to an increase of the polariton splitting. 
Therefore, a weak disorder can increase the robustness 
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FIG. 3: (Color online) Absorption spectrum in the vicinity 
of the lower polariton energy calculated exactly (blue) and 
within the perturbation theory and broadened by Lorentzian 
(red) with 7 = 5fieV. 40 Realizations of energy disorder 
with variation a u = 1 meV on the lattice with 20 2 sites were 
considered. 



and visibility of the polariton splitting with respect to 
the ideal case. The fact that a weak exciton disorder in- 
creases the observability of a coherent effect may be seen 
as counter-intuitive, but it has been also found in similar 
contexts, e.g. in the exciton Aharonov-Bohm effect.— ^ 
Additionally, by increasing disorder the UP and LP peaks 
broaden. Within the perturbative approach we find that 
the broadening of the polariton modes is equal to in 

the lowest order (see Eq. ([29])). the second order being 

proportional to (see Eq. (|3"5"j) ). Moreover, in order to 
check our perturbative approach we have also performed 
a direct comparison between the full numerical solution 
and the analytical result from perturbation theory for 
weak disorder, shown in Fig. [3] and an excellent agree- 
ment is found. 

Notice that for both T and M polaritons, a new cen- 
tral very broad excitonic peak appears. However, in the 
M-point case a central peak corresponding to a pure pho- 
tonic mode, clearly insensitive to the disorder, is super- 
posed to the broad central excitonic background. As the 
width of the energy distribution due to the inhomogc- 
neous energy disorder approaches the value of the lower 
and upper polariton splitting, the three peaks merge and 
the spectrum consists of a central narrow (photon-like) 
peak and a broad excitonic background. 

From an experimental point of view, state-of-the-art 
quantum dot samples have inhomogeneous broadening 
of several meV, which is comparable to typical values of 
the polariton splitting. This makes the observation of the 
polariton splitting challenging. On the other hand, impu- 
rities have very narrow energy distribution and the pla- 
nar lattice can be realized by patterned ion-implantation 
methods. Array of impurties with deep excitonic lev- 
els could be obtained by this method. Assuming deep 
energy levels, Frenkel-like excitons bound to impurities 
will have very small inhomogeneous energy broadening. 
Nevertheless, using ion implantation, it is hard to con- 
trol the number of impurities on every lattice site. Thus, 
the lattice of impurities has a very narrow on-site energy 
distribution but large oscillator strength disorder due to 
impurity number fluctuation. The same kind of disor- 
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FIG. 4: (Color online) Polariton absorption spectrum in log- 
arithmic scale for oscillator strength disorder with variance 
from (black) to 50% (magenta) of the average for polaritons 
at two high symmetry points of the first Brillouin Zone: M 
(a) and V (b). Other parameters as in Fig. [2] 



der may occur in atoms trapped in optical lattices. In 
order to investigate this effect, we consider a Gaussian 
distribution of the exciton oscillator strength in the cou- 
pling constants gj q . A Poisson distribution is another 
possibility but we have not investigated that case. The 
absorption spectra for different oscillator strength disor- 
der are shown in Fig. 2] The behaviour of the LP and UP 
peaks for both V and M-point is similar to the case of the 
energy disorder. In particular, we observe a decrease (in- 
crease) of LP (UP) energy, which is again consistent with 
the perturbation theory as seen in Eqs. (|27|) and (|36|) . 
We also remark that even for very large fluctuations of 
the oscillator strength both the broadening and energy 
shift remain smaller with repect to the inhomogencous 
energy disorder. 

Finally, we have also investigated absorption spectra 
in the presence of positional disorder, i.e. deviations in 
the exciton localization sites from the ideal lattice, as 
shown in Fig. [5] Clearly, this disorder plays a role only 
for M-point polaritons, since the T-point does not probe 
the lattice symmetry (Fig. [SJa) . In contrast to the oscil- 
lator strength disorder and to the energy disorder, posi- 
tional disorder leads to a blueshift (redshift) of the lower 
(upper) polariton peak, which implies a reduction of the 
polariton splitting. This is related to the fact that the 
first order perturbation theory is non-zero in this case, as 
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FIG. 5: (Color online) Polariton absorption spectrum in loga- 
rithmic scale for position disorder with variance from (black) 
to 10% (red), 30% (green), and 50% (magenta) of the lattice 
constant and for polaritons at two high symmetry points of 
the first Brillouin Zone: M (a) and F (b). Other parameters 
as in Fig. [2] 

seen in Eq. (|28p . The behaviour of the central polariton 
(CP) peak is more complex. First, it splits into three 
lines for weak disorder. Then the remaining central peak 
splits again into a doublet for stronger disorder. The 
first splitting can be understood easily taking into ac- 
count CP's three-fold degeneracy^ at qM which is lifted 
in presence of disorder and consequently one energy shifts 
upwards (downwards) and one stays in the middle. With 
increasing strength of disorder, the state with q^j mixes 
with neighboring q and doublet appears. 

Moreover, assuming zero exciton transfer, no anyliti- 
cal expressions can be derived because the CP is TV-fold 
degenerate (in the BZ). On the other hand, assuming a 
finite exciton transfer rate tx ^ 0, the TV-fold degener- 
acy decreases to a three- fold degeneracy at the M-point. 
Since these levels are purely photonic, the first order ma- 
trix elements are zero and second order perturbation the- 
ory has to be used. The matrix elements of the second 
order effective Hamiltonian read 

q»w - ,,°(p)y K'fc'l 2 ,oo^ 

A CPk M - A n'k' 

where 

o(v) \ o(p) *o(r>) * /on\ 

KrZ' = Z^^k+Kk'Vk+K'k'Vmk+KVnk+K'. (39) 

KK> 
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The correction to the unperturbed energy is found by 
averaging the eigenvalues of O, whichs cannot be done 
analytically. A good approximation consists in averaging 
first the matrix elements 



Mp)\ _ ai°(p) 

T"mn I — u kk'KK 



k+K v nk+K' 



(40) 



before calculating the eigenvalues. We have shown that 
this approximation introduces errors of the fourth order 
in the disorder potential and gives good agreement with 
the exact numerical calculation. 

If the positional disorder is of the same order of the 
lattice constant, i.e. for ctr ~ 0.5a, the lattice becomes 
close to a set of random sites. In this limit, instead of 
averaging the eigenvalues, the random phase e tqRj in the 
coupling constants can be averaged first, which leads to a 
position independent value (ge tqRj ) = g. This results in 
a T-point like spectrum with an upper and lower polari- 
ton as clearly shown for the strongest disorder in Fig. [5^ 
(magenta). 



IV. CONCLUSIONS 



or impurities embedded in a planar microcavity. We have 
focused on Bragg polaritons occurring at the Brillouin 
zone boundary and we have shown that small energy dis- 
order, i.e. diagonal disorder, with respect to the light- 
matter coupling is necessary to observe polariton features 
in the absorption spectra. Moreover, we have found that 
polariton is very robust against oscillator strength and 
position disorder. Our results show that a lattice con- 
sisting of deep impurities with Frenkcl excitons or atoms 
in an optical lattice structure are promising candidates 
for the experimental observation of the Bragg polaritons. 



Acknowledgments 



We have investigated numerically and analytically the 
role of disorder in structures consisting of a lattice of dots 



This research was supported by the National Science 
Foundation, Grant No. DMR-0608501. 



1 J. Kasprzak, M. Richard, S. Kundermann, A. Baas, 
P. Jeambrun, J. M. J. Keeling, F. M. Marchetti, M. H. 
Szymanska, R. Andre, J. L. Staehli, et al., Nature 443, 
409 (2006). 

2 R. Balili, V. Hartwell, D. Snoke, L. Pfeiffer, and K. West, 
Science 316, 1007 (2007). 

3 C. W. Lai, N. Y. Kim, S. Utsunomiya, G. Roumpos, 
H. Deng, M. D. Fraser, T. Byrnes, P. Recher, N. Kumada, 
T. Fujisawa, et al., Nature 450, 529 (2007). 

4 H. Deng, G. Weihs, D. Snoke, J. Bloch, and Y. Yamamoto, 
PNAS 100, 15318 (2003). 

5 V. Savona, C. Piermarocchi, A. Quattropani, P. Schwendi- 
mann, and F. Tassone, Phase Transitions 68, 169 (1999). 

6 J. Keeling, F. M. Marchetti, M. H. Szymanska, and P. B. 
Littlewood, Semic. Sci. Tech. 22, Rl (2007). 

7 E. Peter, P. Senellart, D. Martrou, A. Lemaitre, J. Hours, 
J. M. Gerard, and J. Bloch, Phys. Rev. Lett. 95, 067401 
(2005). 

8 J. P. Reithmaier, G. Sek, A. Loffler, S. K. G Hofmann, 
S. Reitzenstein, L. V. Keldysh, V. D. Kulakovskii, T. L. 
Reinecke, and A. Forchel, Nature 432, 197 (2004). 

9 T. Yoshie, A. Scherer, J. Hendrickson, G. Khitrova, H. M. 
Gibbs, G. Rupper, C. Ell, O. B. Shchekin, and D. G. 
Deppe, Nature 432, 200 (2004). 

E. M. Kessler, M. Grochol, and G Piermarocchi, Phys. 
Rev. B 77, 085306 (2008). 

S. Francoeur, J. F. Klem, and A. Mascarenhas, Phys. Rev. 
Lett. 93, 067403 (2004). 
2 A. Muller, P. Bianucci, C. Piermarocchi, M. Fornari, I. G 
Robin, R. Andre, and C. K. Shih, Phys. Rev. B 73, 
081306(R) (2006). 



T. Shinada, S. Okamoto, T. Kobayashi, and I. Ohdomari, 
Nature 437, 1128 (2005). 

F. J. Ruefi, W. Pok, T. G G. Reusch, M. J. Butcher, 
K. Eng, J. Goh, L. Oberbeck, G. Scappucci, A. R. Hamil- 
ton, and M. Y. Simmons, Small 3, 563 (2007). 
I. Kuljanishvili, C. Kayis, J. F. Harrison, C. Piermarocchi, 
T. A.' Kaplan, S. H. Tessmer, L. N. Pfeiffer, and K. W. 
West, Nature Physics (2008), in press. 
D. Gerace and L. G Andreani, Phys. Rev. B 75, 235325 
(2007). 

Y. D. Chong, D. E. Pritchard, and M. Soljacic, Phys. Rev. 
B 75, 235124 (2007). 

F. Bariani and I. Carusotto, Journal of the European Op- 
tical Society 3, 08005 (2008). 

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature (London) 438, 197 (2005). 
K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 
Firsov, Science 306, 666 (2004). 

G. F. Quinteiro, J. Fernandez-Rossier, and G Piermaroc- 
chi, Phys. Rev. Lett. 97, 097401 (2006). 

R. Houdre, R. P. Stanley, and M. Ilegems, Phys. Rev. A 
53, 2711 (1996). 

V. M. Agranovich and Y. N. Gartstein, Phys. Rev. B 75, 
075302 (2007). 

J. J. Hopfield, Phys. Rev. 112, 1555 (1958). 
G. D. Mahan, Many-particle physics (Kluwer Aca- 
demic/Plenum Publishers, New York, 2000). 
M. Grochol, F. Grosse, and R. Zimmermann, Phys. Rev. 
B 74, 115416 (2006). 



M. Grochol and R. Zimmermann, Phys. Rev. B 76, 195326 
(2007). 



